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Abstract 

We investigate the weakly nonlinear evolution of a triad of nearly-neutral modes su- 
perimposed on a mixing layer with velocity profile u = u m + tanhy. The perturbation 
consists of a plane wave and a pair of oblique waves each inclined at approximately 60 to 
the mean flow direction. Because the evolution occurs on a relatively fast time scale, the 
critical layer dynamics dominate the process and the amplitude evolution of the oblique 
waves is governed by an integro-differential equation. The long-time solution of this equa 
tion predicts very rapid (exponential of an exponential) simplification and we discuss the 
pertinence of this result to vortex pairing phenomena in mixing layers. 
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1. Introduction. 


Transition to turbulence in unbounded, incompressible mixing layers occurs in several 
stages. The initial instability, according to most experiments, is primarily two-dimensional, 
an observation consistent with linear stability theory where the largest growth rates cor- 
respond to plane waves propagating in the mean flow direction. Nonlinearity results in 
eventual saturation of the linearly most amplified wave and as this mode tends toward 
equilibration, a subharmonic wave having one-half the frequency of the initial disturbance 
makes its appearance (as do higher harmonics). This intriguing phenomenon is observed 
in smoke visualization experiments and is also apparent from hot-wire signals. It is usu- 
ally described as “vortex pairing”, i.e., the initial disturbance manifests itself as a row of 
vortices and adjacent pairs of vortices begin to roll around each other, finally coalescing 
into a single vortex. 

The first experimental observations of vortex pairing were reported around 1960 and 
concerned laminar flows undergoing transition. It was not until the mid 1970s, however, 
that interest in the subject greatly expanded because of the first reports that this process 
occurs in turbulent mixing layers, as well, and is a dominant factor contributing to the 
downstream growth of the mixing layer. A comprehensive review of these experiments, as 
well as the early theoretical developments, can be found in the survey article by Ho and 
Huerre (1984). The same article also discusses some early observations of three-dimensional 
structures which take the form of counter-rotating vortices whose axes are in the direction 
of the mean flow. Considerable experimental work has been reported more recently on 
the three-dimensional structures. Huang and Ho (1990), for example, discuss spanwise 
wavelength doubling due to merging of the counter-rotating streamwise vortices. 

The present paper is, of course, concerned with the mathematical analysis of these 
phenomena. A principal contribution toward that end is contained in ideas put forward by 
Kelly (1967) who showed that certain features observed in the early stages of vortex pairing 
could be modelled (and in some cases predicted) by studying the resonant interaction of 
two disturbances whose wavenumbers are in the ratio 2:1. The shorter wave in Kelly’s 
analysis was taken to be periodic in space and time (i.e. neutral) whereas the long wave 
had, initially, a much smaller amplitude. Its subharmonic instability corresponds to the 
onset of pairing. Subsequently, Patnaik et al. (1976) performed numerical simulations of 
instability in a stratified mixing layer and showed that the entire pairing process could 
be modelled by following the temporal evolution of two interacting Fourier modes with 
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wavenumbers related as in Kelly’s analysis. 

Patnaik et oil. discovered by trial and error that the interaction of two such modes was 
greatly influenced by the phase relationship between them and in one limit it is very weak 
(they termed this a “shredding interaction”). This sensitivity to the relative phase is in 
fact a feature of several weakly nonlinear resonant interaction theories and had been noted 
earlier in the context of capillary-gravity waves, as pointed out by Collins and Maslowe 
(1989). There are some potentially important technological consequences of the foregoing 
remarks because they suggest that by forcing the interacting modes with the desired phase 
relationship the rate of growth of the mixing layer can be controlled. This has been verified 
experimentally by Hussain and Hussain (1989). 

The role of the relative phase was emphasized in an extension of Kelly’s analysis by 
Monkewitz (1988) in which the fundamental mode was taken to be nearly neutral, whereas 
the subharmonic was spatially amplified. Even though the analysis is described as weakly 
nonlinear, the wavenumber a of the subharmonic was taken to have an imaginary part as in 
linear spatial stability theory because the author argues that the amplification factor, —<*i, 
is too large to be accounted for by introducing a slow scale in the flow direction. Despite 
this inconsistency, comparison of some experiments with an “ad hoc” version of the theory 
taking into account weakly nonparallel effects shows reasonable qualitative agreement. 

We do not believe that this is entirely fortuitous and in the present paper seek to 
formulate a rational analysis of resonant triads in a tanh y mixing layer. There are two keys 
to overcoming the obstacles cited by Monkewitz. The first is to employ a pair of oblique 
waves as the subharmonic so that all modes are neutral or at least nearly- neutral in an 
asymptotic sense. The resulting analysis which is related to studies by Wu (1992) of the 
Stokes layer and Goldstein and Lee (1992) of the adverse pressure gradient boundary layer 
is greatly facilitated in the example of the tanh y mixing layer by the availability of a closed- 
form neutral solution. Consequently, we are able to determine explicitly the coefficients 
of all terms in our amplitude equations without the need for numerical computation or 
long- wave expansions. 

The second key is to use an appropriate analysis of the critical layer (a thin layer 
centered on the point y c where the phase speed of the interacting waves is equal to the 
mean flow velocity) because the critical layer dynamics seem to dominate the evolution 
of the flow for a significant distance beyond the region of exponential amplification of the 
linearly most unstable wave. This was demonstrated convincingly by Hultgren (1992) who 
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employed a quasi-equilibrium nonlinear critical layer analysis downstream of the initial 
region of linear instability. In addition, weakly nonparallel effects were taken into account 
outside of the critical layer. Excellent agreement was obtained with several experiments 
up to the onset of vortex pairing. Briefly, one reason why such favourable comparisons 
can be achieved despite the questionable small-amplification assumption in the theory is 
that due to the spreading of the mixing layer the wavenumber on a nondimensional basis 
increases as the shear layer thickens and it approaches the neutral value near saturation 
(a n = 1 for a tanh y mixing layer). 

An extension of Hultgren’s analysis to include resonant interactions would in all like- 
lihood lead to the best agreement with the experimented data available to date. However, 
because simulating these experiments is not our primary objective, we will employ a pro- 
cedure that is appropriate to the early stages of pairing and that will facilitate comparison 
with numerical simulations of the Navier-Stokes equations. Specifically, we will employ 
what some authors have termed a non-equilibrium critical layer in which the flow evolves 
on a relatively fast time scale, i.e., e l / 4 t, where e is an amplitude parameter, rather than 
e}! 2 t which would be the slow scale for a nonlinear critical layer (see, e.g., Maslowe (1986) 
for a review of critical layer theories). Although a spatial theory using a slow e l ^ 4 x scale 
would be appropriate for comparison with experiments, most numerical simulations are 
based on the temporal theory for reasons of economy. In any case, as shown in Section 5, 
our amplitude equations can be converted easily to the spatial case. 

Due to the above choice of slow time scale, our critical layer equations can be char- 
acterized as unsteady and weakly nonlinear. According to the amplitude equations that 
result, the nearly-neutral plane wave will amplify exponentially as in linear theory, but 
the oblique mode amplification is governed by an integro-differential equation of the sort 
first derived by Hickernell (1984) in the context of a mixing layer on the beta plane. The 
asymptotic solution of this equation for very large times is shown in Section 4 to exhibit 
exponential of an exponential amplification. Finally, in Section 5, we compare our results 
with some numerical simulations and also discuss their possible relevance to experiments. 

2. Formulation and outer expansion. 

We consider the stability of the dimensionless mixing layer profile 

u(y) = u m + tanh y (2.1) 


4 


by adding a small perturbation of 0(c), where e 1 is a dimensionless amplitude param- 
eter. The equations of motion can be written 

|3 + (q-V)q=-Vp (2.2) 

and 

V • q = 0, (2-3) 

where we have supposed the fluid to be inviscid and incompressible. The temporal stability 
problem is independent of u m , the mean velocity in (2.1), so we will set u m = 0 until later 
in the paper when we discuss the spatial case. 

The velocity components q = (u + eu,ev,ew) and the perturbation pressure ep are 

expanded as follows: 

u = u< 1 >+£ 1 / V j >+£ 1 'V 3 > + --- (2-4) 

t, = »<>> + £ 3 / V 2 > + £ 3 / V 3 > + •• ■ (2.5) 

w = w f,) + + e'/V 3 ’ + ■ • ■ (2.6) 

and 

p = p( 1) + £ l/y 2) + £ l/2 p (3) + . . . . (2.7) 

In the linearized problem, perturbations proportional to exp[i(ax+/?z— act)] are considered 
and the y-dependent part of tq, which we denote tq, satisfies the Rayleigh equation 

(u - c)(tJ" - a 2 v x ) - u "tq = 0, (2.8) 

where a = (a 2 + 0 2 ) 1 / 2 . 

For the mixing layer profile u — tanhy, Curie (1956) has found the following neutral 
solution of the eigenvalue problem consisting of (2.8) and the boundary conditions that 
tq — > 0 as y — ► ±oo: 

tq = sech y, a 2 = a 2 + /? 2 = 1 and c = 0. (2-9) 

We employ this basic solution to construct an initial disturbance of the form 

u (1) ={A 20 (T)e iai + A^e-'^Mv) 

+ {i n (T)e ,ai/2 + A\ x t- iax/2 } tq(y) 2 cos 0z, (2.10) 
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where v\ is given by (2.9) and T is a slow time scale. The disturbance can be described 
as a plane wave plus a subharmonic consisting of two oblique waves of equal amplitude 
and equally inclined with respect to the mean flow so that in the spanwise direction we 
have a standing wave. The factor of 2 is introduced because it will simplify the subsequent 
development to utilize complex exponentials along with the identity 2 cos f3z = exp (i/3z) -f 
exp (—i0z). 

A triad of neutral modes satisfying the resonance conditions exactly could be chosen 
according to (2.9) if a = 1 for the plane wave and the subharmonic is comprised of the 
oblique modes having a = 1/2 and /? = ±\/3/2. We depart only slightly from this scheme 
by allowing the plane wave to be weakly amplified and set 

a — 1 — e l t A a\ . ( 2 . 11 ) 


From Lin’s perturbation formula, the linear growth rate is 

aci = — (1 — a) = — e 1 / A ai 
7 r 7T 

and it follows that an appropriate slow time scale is T = e 1 / 4 t. 

The procedure for determining the remaining terms in (2.4) - (2.7) will be outlined 
only briefly because the details are straightforward and in many respects parallel the in- 
vestigation of Benney (1961). (A principal difference is that Benney employed a steady, 
viscous critical layer. This, however, has little effect on the first terms of the outer expan- 
sion.) When a pair of oblique waves is included in the basic disturbance, then both 
and u/ 1 ) contain terms having first-order poles at the critical point y = 0. For example, 

w yi = cosech y, (2-12) 


(cf. Section 4 of Benney) where we introduce the notation w ^ to mean the term at 
0(£^” -1 ^ 4 ) containing the factor exp{i(^x + 'J%mz)j2}. 

The presence of singular terms such as that in (2.12) generates discontinuities in 
certain higher-order velocity components and these must be smoothed out by the critical 
layer solution. Although the general form of the amplitude equation is determined by the 
critical layer, the values of coefficients appearing in those equations must be determined 
by matching to the outer expansion. 

Returning now to the expansions (2.4) - (2.7), at 0(e}t A ) the quantity satisfies a 
nonhomogeneous Rayleigh equation which can be written 


C\v = — 2aiAnsech y — 4iA , 11 sech 2 y cosech y, 


(2.13) 
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where C x is the Rayleigh operator in (2.8) with 3=1. Because both homogeneous solutions 
axe known in closed form, the method of variation of parameters can be used to write the 
general solution of (2.13) as 


v n = Cji }± sech y + D[\ }± (y sech y + sinhy) - <*1 A u coshy 


(2.14) 


+ 2iA\ 
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2sech y / yicosech 2yi dy\ — (y sech y + sinhy) log(tanh |y|) , 

Jo 


where ± denotes above and below the critical layer. Imposing the homogeneous boundary 
conditions as y — > ±00 shows that there is a jump in which is related to the amplitude 

of the oblique waves by 

D [\ | )+ - D { ^~ = 2a!A u . (2.15) 

/rt\ /2\ (2) 

Equations (2.2) and (2.3) can now be used to solve for tr n ; , and p u . 

We next repeat this procedure for the 0(e 1/4 ) terms involving the plane wave. The 
quantity satisfies 


C \ = — 2 oiA 2 osech y — 2 iA 2 0 sech 2 y cosech y, 


(2.16) 


where C\ again is the Rayleigh operator with a = 1. Using variation of parameters, we 
obtain the same solution as in (2.14) except that the factor 2 in front of the term in 
square brackets is absent. Imposing the boundary conditions as y — *• ±00 yields the jump 
condition 

D ( 2 2 o H ~ £>2 2 o )_ = Sc^o. (2-17) 

It will be seen that matching the jumps in (2.15) and (2.17) to the critical layer 
solution derived in the following section will lead to the amplitude equations governing the 
temporal evolution of A\\ and A 20 . Before continuing, however, we will note one procedural 
difference compared with Hickernell (1984) and other papers employing non-equilibrium 
critical layers. Because an analytical solution was not available for the outer problem, the 
jump conditions in these papers were derived from a generalization of the usual adjoint 
orthogonality condition in which it is implicit that an outer solution can be found satisfying 
the boundary conditions. Here, on the other hand, we employ the procedure described in 
Section 3 of Benney and Maslowe (1975) and actually find the outer solution. The two 
procedures should lead to the same result for the amplitude equations. 
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3. Critical layer analysis. 

To obtain evolution equations for A 2 o and .An, we shall now pose an inner expansion 
in the critical layer in order to obtain expressions for the jumps. To this end, we introduce 
the inner variables Y = £ -1 / 4 y, U = V = e -1 ^ 2 u, W = e -1 / 4 io and P = £ -1 p 

inside the critical layer, where it should be recalled that e was the order of magnitude of 
the disturbance in the outer expansion, so that the governing equations become 


U T + uu x + VU Y + wu z + e~ 1/2 P x = 0 (3.1) 

V T + UV X + VV Y + WV Z + Py = 0 (3.2) 

W T + UW X + VW Y + ww z + e~ 1/2 P z = 0 (3.3) 

U X + V Y + W Z = 0. (3.4) 


The form of the outer solution written in the inner variables suggests that the inner 


expansion is of the form 

U = Y + s l ' 2 U 2 + £ 3 / 4 Z7 3 + eUi + ■ ■ • (3.5) 

V = £ 1/2 v 2 +£ 3/4 V r 3 +eV 4 + --- (3.6) 

W = e l/2 W 2 +e :ili W 3 +eW4+--- (3.7) 

P = P 0 + £ 1/4l Pi + e 1,2 P 2 + e 3 / 4 P 3 + eP 4 + • • • . (3-8) 

Substituting this expansion into the governing equations, collecting powers of £ and 
grouping terms with the same x dependence, we arrive at a series of equations of the form 

nY 

4'T + i^-4’ = x(Y,T) (3.9) 

which have solutions of the form 

4>= f T x(^To)e n ^( To - T >dT 0 . (3.10) 

J — oO 


We will also need to calculate the jumps in several quantities across the critical layer with, 

f°° 

for example, the jump in <p given by / when we come to do this, it should be 

J — OO 

recalled from the definition of the Fourier transform that, for real a^O, 

/ e iaYT dY = — 6(T). 

J - 00 a I 
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For the inner expansion, we shall use the notation Um,n to mean the term at O (e^ 4 ) 
accompanying exp (io {rnx + fiy/Sz] /2); in what follows, we will only give the terms cor- 
responding to non-negative m and n, with the remaining terms following from symmetry, 
for example, U ^2 o = ^ 20 ^ • Additionally, only those terms necessary to evaluate the jumps 
will be determined. 


3.1 0(e 1 ^ 2 ) terms. 

At this order, we find that, in addition to both the two-dimensional and the oblique waves, 
there is a non- zero component of the mean flow present in the x-direction, Uqq = —Y / 3, 

with v 0 ( 0 2) = w<? = 0 . 

For the two-dimensional wave, we find that U-J, 1 — W 2 o^ = 0 and V 2 o = ^ 20 , with 
the relevant pressure terms given by 


P 2 ( o 0) = iM 0 

p(l) *>o(2)(i) 

-*20 ” lo 20 


(3.11) 

(3.12) 


p(2) _ 
*20 — 


iA,„Y 2 


- A 20 Y + S.-4. + + iCg** - 2 (3.13) 

Since P^ and must be continuous, this tells us that 
and do )(+) + 2z^ 2 0 )(+) ' = C^ )(_) + 2i D^~ Y . 


'20 

For the oblique waves, we find that Vn = An and that satisfies 


V 


( 2 ) 

11, T 


, iY M 2 )_ 3 A 

+ Y Ul1 ~~4 Al1 


with solution u [ ^ = — — f An(To)e 2 ^ To T ^dTo 

^ J — oo 

and W™ = = y T f T A„(T„)e i fm-l’)dr, 

** J — oo 


V5 

with the relevant pressure terms given by 


(o) _ i An 
*11 — 


p(l) 

*11 “ 


iC\\ 


< 2 )(±) 


(3.14) 

(3.15) 

(3.16) 

(3.17) 


and = 


i? = --^1 - A n Y + MAu + 


l OL I C 


r(2)(±) 


iC 


(3)(±) 


+ 2 D[\ )(±) '. (3.18) 

4 2 2 11 v 

Since Pfp and P x (2) must be continuous, this tells us that Cj 2)(+) = ^ = C ^ and 

c!i )(+) + 4iD^ K+) ' = CSJ )(_) + 4iD[] )( ~ y . 


9 



3.2 O (s 3 / 4 ) terms. 

At this order, only the exp(mx) terms and the exp (ia (x + y/Z z) /2) terms need to be 
calculated. For the two-dimensional wave, we find that C/ 20 ' = = 0 and that V 2 o == 

C$f -cnA 20 (T). 

For the oblique waves, we find that ; = Cf/ — aiAu/4 and that 


w( 3 ) i *— 
w n t T + 2 rr 11 


w d) = vgMn j VScff iajYU™ 


+ 


4 


2y/Z 


VzclV Vz ai r T x' 11 (r, ) «^( r *- T )<ni 


£ 


which has a solution 




( 3 ) 


= : T f 1 C;; ) (To)e^( T °- T >dTo 

^ J — OO 


11 


+ 


V^CKi 


f T (T 0 - T)A; i (To)e i f( T ‘»- T )dTo, 
J — OC 


and that 


jr/( 3 ) i —Tj^> — 

u n,YT t 9 l/ u,y — 


»r 


( 3 ) _ 


iVSHtf’ ictYV^r iai U <? 


+ 


which has a solution 


= | f C^(To)e^< r °- T >dT 0 
- ^ f T r A' 11 (To)e^( r °- r )dTodr 1 

^ */ — oo %/ — OO 

f C’i? ) (T 0 )e^( T °- T >dTo 

- ^ f T (To -T)A' n (To)e^( To - T >dTo 

** — OO 


(3.19) 


(3.20) 


(3.21) 


(3.22) 


so that ufi* + y/ZWff = 0. 


3.3 O (e) terms. 

It is at this order that we shall evaluate the jumps necessary to obtain the evolution 
equations. For the two-dimensional wave, from the outer, we know that as Y — * ±oo, 

t/ 2 (4) + iA 20 {T)Y ~ 2A 20 log | y | +A; 0 + 2 iZ>g )(±) + O (F _1 ) 
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so that 


-yjy^y 


(3.23) 


We find that satisfies 


(j^ + *^ U^Y = ~ iP 20 ,Y + ia l P 20 ,Y - ^ iU YL U YL,Y ~ 2 A U U 1 i]y,Y 


= -- F F /" Tl A 11 (To)A 11 (T 2 )e i r (To+T2 " 2T) dTodT 1 dT 2 

8 J — oo J — oo J — oo 


-ILLL 

3 ,r rTj ,r, 

^ J— CO J — oo J — oo 


^n(TMii(To)e^ (r °- T) dTodT 1 dT 2 


+ ^ 20 (T)y + iA' 20 (T) 


(3.24) 


so that 
^ 2 (>)y *-^20 

o />T /"Ti 


=- F F (-3T 2 + 3TTo - T 0 2 + 3TTi - T 0 Ti - T^Ai^ToMi^TOe ’^ (To+Tl 2T) dT 0 dT i 

8 7-00 7-00 

+2i F A 20 (To)e iYiT °- T) dTo 

J — OO 

= -| r° /"(T 0 2 + 3 Tor 1 + 3 r 2 )A 11 (T 0 )A 11 (T 1 )e- i f (ro+ 2 ri) <iTodT 1 
8 7o 7o 

+2i r A 20 (T - T 0 )e- iYr Mr 0 . (3-25) 

Jo 

The jump in (^^(V', T) + 2^20 across the critical layer tells us that 

1 roc v 

^ = jL (°& + ^ 20 ) ' 

= — F° l°° /°°(r 0 2 + 3r 0 ra + 3r 2 ).4 u (T - r 0 - r 1 )A u (T - n)e- ^ T o+2ti) dr 0 dT X dY 

32a x 7 -oo Jo Jo 


32a n 7 - 


-1 /»oo roc 

+ 7-/ / A M (T-ro)e- iYr ‘dr 0 dY 

2c*i 7-oo 7o 


= _^!!L r / 00 ^(To+2r 1 )(r 0 2 + 3ror 1 +3r 1 2 )Aii(T-To-r 1 )^ 11 (T-r 1 )drodr 1 

16 ^i 7 o 7 o 

rOO 

+ — / 6 (r 0 ) A 20 (T - T 0 )dT 0 

<*1 7 o 


= ^^o(T), 
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and hence one equation is 

A 20 — — A 2 0 = 0. (3.26) 

Za 1 

Likewise, for the oblique waves, from the outer, we know that as Y — * ± 00 , 

U[\ ] + V3 W}? + 2iA n Y ~ 8A' U log I F I +4 A u + 4i£>< 2 1 )(±) + O (F -1 ) , 


so that 


in 


= 2^( D 


( 2 )(+) 


D 


Wi- 
ll 


’) 


= sk f-„ + ^'■Y + 2iAi1 ) iY ■ 


We find that 


so that 


( 


Yt + y) + ^ W wy) 

= -2iP$y + 2 iajPtfy - 2A 20 (T)u} 2 1 >yy - A n (T)U$yy 


= YA n (T) + 2iA n (T) 

rT r T 2 i*Ti 


- | / f 2 [ 1 ^2o(T)AI 1 (To)e^(- To+T >dTodT 1( fT 2 

^ J — OO J — OO J — OO 


+ 2* A u 


ii,y 
rT 


(3.27) 


=4 i [ A' u (T 0 )e‘f< T °- T >dTo 
J — OO 

- i / T r ( T ° - ri) 2 A 20 (T 1 K 1 (To)e i ?(- T ‘> +2r ‘- T )dTodT 1 
® v/ — 00 J — OO 

roo # ;y> 

= 4 i A U (T - r 0 )e~ dr 0 

Jo 

o y*oo roo 

_ 8 i 0 L r ° 2A2o(T " W - Ti - ro)e i ?( ro -^)drodr 1 . (3.28) 

The jump in ^£7ii^(F, T) + (Y,T) + 2iAnF^ across the critical layer tells us that 
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' 4u = st; /” ( a “’ r + + 2iAu ) dY 

= — r r A ^ T - r ^ e-^drodY (3.29) 

2ori J_o o ./o 

+ TT- [°° r / 00 ^^2o(T-r 1 Mt 1 (T-r 1 -r 0 )e^ ro - ri) rfr o( ir 1 dr 
04a! y 0 ./o 

O /-oo 

= — / <5(ro)A u (T-r 0 )(fro 
«x 7o 


+ 


3i7r 

16a 


n oo 

ro 2 

. 


6 (n - t 0 ) A 20 (T - r^AhiT - t x - To)dTodr l 


= -A n (T) + r^AMT - to )A’ U (T - 2r 0 )dr 0 , 
ai 16aj J 0 

and hence the amplitude equation for the oblique waves is 

o ■ roo 

A' M --An=-fi / rlA 10 (T-r 0 )A\ l (T-2T 0 )dr a . 

7T 10 Jo 


(3.30) 


4. Solution of the amplitude equations. 

Equation (3.26) has the obvious solution 

A 20 (T) = a 20 e 2aiT/7r , (4.1) 

where a 2 0 is a constant, and substituting this result into (3.30), we obtain 

A„(T)-— A' n (T)=^^ /“t 0 2 « 2 «'< t -'*WM; i (T-2to)* o1 (4.2) 

a i ibai 

the left-hand side of this equation suggests that we write -4n(T) = Bu(T)e aiT /* , which 
leads to an equation for B n (T ) which can be cast in the form 

B'u = r (T-rof^'^’BUMdro, (4.3) 

120 7- oo 

which we can integrate with respect to T to get 

Bu(T) = “u-‘^ f ( T - r„) 3 (4.4) 

where a\\ is a constant. 
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4.1 Solution by iteration. 

A standard technique for solving equations such as (4.4) is the method of successive ap- 
proximations (Hildebrand, 1965, §3.9), whereby we solve the equation iteratively, taking 
= an to be our initial guess, and then for successive iterations talcing 

rp 

= <„ - ^ jf (T - r 0 ) 3 (4.5) 

which leads to a solution in the form of a series 


Bu(T)= 


(4.6) 


m= 0 


with /3 0 = <Xn and /3 m = — 3ia2o/3^ n _ 1 Tr i /(2 10 m i a\). The ratio test tells us that this series 
converges absolutely for all T < oo. We can write this series for Bn as 

Bn(T) = i (an + *;,.<■*•->) f ^ 

+ I (a,, - >) E (4-7) 


It follows that in order to evaluate Bn, we must evaluate the series E x m /(m!) 4 . We 
know (Abramowitz Stegun, 1964) that 


m=0 


.2 / s ^ (-ir(2m)!(^ 2 r 

7 ° W= E idi ! 


m=0 


(m!) 4 


and 


(2^j! = s/ c ( - pr2 "' le “ r<ip ’ (48) 

the latter of which is Hankel’s contour integral, where C starts at +oo on the real axis, 
circles the origin in the counterclockwise direction, and returns to the starting point. It 
follows that 


oo (l~2\m oo 

E Wf l sr 

( m !) 4 ^ 

m=0 v J m — 0 

so that 


Qx 2 ) m (2m )! 
( m !) 4 


X 




»„m - - 4 («. . jf a . 


ia^e 


(4io) 
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4.2 Asymptotic solution for large T. 

In order to evaluate the large T behavior of (4.10), it is necessary to evaluate 

Lir J KiJ dp 

and , . 

IcT II (l) dp 

as x — * oo. If we set p = qx in (4.11), it becomes 


Since w(x,q ) = -qx - log q has a saddle point at q = -1/x, we can evaluate the large x 
behavior of (4.13) by setting 


dp 

(4.11) 

dp 

(4.12) 


(4.13) 


1 iqo 
^ x x 


which tells us that 


j - — J e ? ° /2 J$ (~ x ) dqo ~ iey/^Jl ( x ) ■ 

Following a similar procedure with (4.12) tells us that 

J ——— Iq f ^ dp ~ ieV^Tr/o (x) . 

Now, as x — *• oo, we know (Abramowitz &: Stegun, 1964) that 

Jo(x) ~ — (sin2x + 1) 

7TX 


tf(*) 


so that 


f tZji (i) dp „ + 

Jc p \pj V* x 

I L z n (1) dp „ hfipl, 

Jc p \pj v* x 

and therefore our expression for Bn(T), namely (4.10), tells us that as T — ► oo, 

2 5 /Ve(a„ +ia; 1 e i “«<“">) «, T/ » (, , ... Wl a 2 o 


_ 2 s / 2 tt;e(a 1 i+. n ; 1 e i »'«l«“)) :n , T; , C 
7T 7 / 2 a/ 3 ^/l °20 | \ 

2 5 / 2 a?e(a 11 -^at 1 e 1 ^ tt ^) ;a1^/ , 

7T 7 / 2 \/3a/| a 2 o I 


+ sin 


^ 8a? 

n/ 3^/1 a 2 o [^ 2 e air/7 
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Since we set An(T) = B\\(T)e aiT ''* , we have an expression for An as T — ► oo, 


^u(T) 


2 5 * / 2 afe(an + ia^e* ar 6( a 2°)) ( ^ . / >/3\/l a 20 K 2 e“ lT / ,r 

~ — ” — , — = I 1 -j- sin I — — 0 

7T 7 / 2 \/3\/| 020 | \ \ 8a l 

2 5 / 2 a 2 e(a n - ia^e^ 1 *^) / V^a/I a 20 |7r 2 e ai:r / T 

+ exp ( R 


(4.20) 


If the coefficient in front of the second term in (4.20) is non- zero, an ^ then 

the second term will dominate the large T behavior, and we have 

, 2 5 / 2 a 2 e(a 11 -iat 1 e iai 8(“^) ( V3^/\^~\n 2 e^ T /” 

All(T) — exp ( si 

and we will have so-called e &T growth. However, when an = ia^j e' ar s( a 2 °) , so that the 
second term in (4.20) vanishes, we are left with 



An(T) 


2 7 / 2 a\t 


7tV 2 v/VI °20 


1 + sin 


V*V\ a 20 | n 2 e aiT / 

8a 2 



(4.22) 


so that as T — * oo, the oblique waves stop growing and rather oscillate about a fixed value. 
This corresponds to the shredding interaction discussed in Section 1. 


5. Concluding remarks. 

In the preceding sections, we have investigated the weakly nonlinear interaction of a 
triad of modes comprised of a plane fundamental and a pair of oblique subharmonics that 
in the linear neutral limit satisfy the resonance conditions exactly. This has enabled us 
to carry out a rationed asymptotic analysis of the interaction; when only two-dimensional 
waves axe included this cannot be done. Singularities that occur at the critical point have 
been resolved using a non-equilibrium critical layer that evolves slowly because the plane 
wave is weakly amplified on a linear basis. To the order of our analysis, the exponential 
amplification of the plane wave is not altered by the interaction and, for that reason, we 
term this a parametric resonance. The common amplitude of the oblique waves, on the 
other hand, satisfies an integro- differential equation, namely, (3.30). This is in contrast to 
the more conventional resonant interaction theory which employs a viscous critical layer, 
an et slow time scale instead of the present and leads to a pair of coupled nonlinear 

ordinary differential equations [see, e.g., eqs. (1) of Collins and Maslowe (1988)]. 
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Both an exact and a long-time asymptotic solution of (3.30) were found in Section 
4. Of particular interest is the asymptotic result (4.21) showing very rapid growth of the 
oblique waves for T 1. This raises the question of the extent to which such amplification 
might be observed in a numerical solution of the full Navier-Stokes equations. Although 
many numerical simulations of unstable mixing layers have been reported in recent years, 
none employs initial conditions that would permit the resonance studied here to be ob- 
served. However, several computational papers are at least indirectly related to our own 
work in that both oblique and plane wave perturbations were superimposed on a plane 
mixing layer. These include articles by Metcalfe et al. (1987), Moser and Rogers (1993) 
and Klaassen and Peltier (1989). In each of these papers, significant three-dimensional 
instabilities were reported, the third differing from the first two in that a Floquet analysis 
was employed. Unfortunately, these papers all view the three-dimensional perturbations 
as secondary instabilities, whereas in our analysis the amplitude of the oblique waves is 
the same order of magnitude as the plane wave. 

Nonetheless, Metcalfe et al. have considered a set of initial conditions that is, in some 
sense, related to our own. In particular, results depicted in Figures 12 and 13 of their paper 
correspond to the initial superposition of a plane wave with a = 0.4, a subharmonic with 
a = 0.2 and a pair of oblique waves having (a, /3) = (0.4, ±0.2). The Reynolds number 
in these computations was 400 which form a stability standpoint is practically inviscid. 
Although the wavenumbers seem, at first glance, to correspond to modes that are too 
highly amplified for comparison with our theory, the thickness of the mixing layer doubles 
during the computation so that at some point they are not far from the linear neutral 
values. 

On the basis of weak nonlinear interactions, it could be argued that the subharmonic 
and oblique modes interact to produce a subharmonic resembling our own except that the 
angle of inclination to the flow direction would be 45° instead of 60°. This discrepancy is 
not serious and could be accounted for by modifying the coefficient of the linear An term 
in (3.30). What really precludes any direct comparison with the present theory is that 
the amplitude of the plane wave in the simulations of Metcalfe et al. is several orders of 
magnitude larger than the oblique and subharmonic perturbations. Even so, after a period 
of time the amplitudes of both subharmonic and oblique modes become comparable to the 
fundamental. This is especially true in Figure 12 which pertains to a case where the oblique 
modes have a much larger initial amplitude than the subharmonic. 
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Given that even if it is not the case initially the amplitude of the oblique waves can 
become larger than that of the plane wave according to our analysis (and the numerical 
simulations do not dispute this) a different scaling in which the oblique waves have a larger 
amplitude may be appropriate. Such a distinguished limit has been identified by Goldstein 
and Lee (1992) in their study of the adverse pressure gradient boundary layer [see, also, 
Wu (1992) for a subsequent application to the Stokes layer]. This scaling leads to coupled 
nonlinear integro-differential equations for the amplitudes which develop a singularity that 
signals the onset of a possibly fully nonlinear stage. The singularity appears, in fact, even 
in the earlier work of Goldstein and Choi (1989) where only the pair of oblique waves was 
present without a plane wave. 

We have formulated this sort of theory for the mixing layer and are presently preparing 
a paper on that topic. The amplitudes of the oblique and plane waves are taken to be, 
respectively, 0(e) and 0(e 4 / 3 ) and the slow time scale becomes e 1 / 3 f instead of the present 
s l / 4 t. These scales are identical to those employed by Wu (1992). When the amplitude of 
the oblique waves is taken to be small compared with the plane wave the new amplitude 
equations reduce to those of the parametric resonance stage, as was shown by both Wu 
(1992) and Goldstein and Lee (1992). 

Finally, some comments about the spatial case are appropriate due to its relevance 
to experiments. As discussed in the Appendix of Kelly (1967), it is necessary that u m be 
greater than zero to have spatial growth and so we discuss the case u m = 1. In a frame of 
reference moving with the real phase speed (if it were the same for both the oblique and 
plane waves), the evolution equations (3.26) and (3.30) would be unchanged except that 
the independent variable would be X = e l ^x instead of T. 

An important difference, however, between the spatial and temporal cases is that 
spatial modes are dispersive for the mixing layer profile (2.1), whereas unstable temporal 
modes always have c r = u m . Numerical computations for the spatial case with u m = 1 
[see, e.g., Figure 3 of Maslowe and Kelly (1971)] show that c r does not differ greatly from 
1 except for relatively long waves with a < 0.4 and this is why numerical simulations of 
the temporal case seem to yield results that can be compared with experiments. In some 
of the experiments reported to date, the spanwise wavenumber (3 is small (although it can 
vary with downstream distance), so the angle of the oblique waves’ inclination is not close 
to the 60° in our analysis. However, provided that the dispersion is not large, the present 
analysis could easily be modified to take into account detuning and this should be done if 
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comparisons were to be made with a particular experiment. 

For the nearly neutral modes considered in this paper, the difference in phase speeds 
in the case u m = 1 turns out to be negligible and could be completely eliminated by having 
the inclination of the oblique waves increased slightly. If we write 

c = 1 — e x t A c \ , (5.1) 

then for the plane wave c\ = Aai/n 2 , whereas for the oblique waves c\ — 4.5c*i/7r 2 . 
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